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The ferrimagnetic spinel C0V2O4 has been a topic of intense recent interest, both as a frustrated 
insulator with unquenched orbital degeneracy and as a near-itinerant magnet which can be driven 
metallic with moderate applied pressure. Here, we report on our recent neutron diffraction and 
inelastic scattering measurements on powders with minimal cation site disorder. Our main new 
result is the identification of a weak (— ~ 1(U 4 ), first order structural phase transition at T* = 90 
K, the same temperature where spin canting was seen in recent single crystal measurements. This 
transition is characterized by a short-range distortion of oxygen octahedral positions, and inelastic 
data further establish a weak A ~ 1.25meU spin gap at low temperature. Together, these findings 
provide strong support for the local orbital picture and the existence of an orbital glass state at 
temperatures below T*. 

PACS numbers: 61.05.fm,75.25.-j,75.25.Dk, 78.70.Nx 


I. INTRODUCTION 

The spinel vanadates, AV2O4, are interesting materi¬ 
als which lay at the heart of several distinct topics in 
condensed matter physics. In recent years, the focus has 
been geometric frustration, and the role played by the 
vanadium orbital angular momentum degrees of freedom 
in relieving the significant spin frustration inherent to the 
spinel structure 1 . In insulating vanadates containing di¬ 
valent A-site cations, the octahedrally coordinated V 3+ 
comprise a pyrochlore sublattice with S = 1 spins and 
an orbital triplet degeneracy in the ideal cubic crystal 
field environment. When the A 2+ cation is non-magnetic, 
this degeneracy is relieved at a Jahn-Teller structural 
transition from cubic-to-tetragonal global symmetry, fol¬ 
lowed by an antiferromagnetic spin order transitionjjpon 
cooling 2 10 . In materials MnV^Cb 11 1§! and FeV204 18 ^} 
the A-site also contains a spin, and A — B superexchange 
interactions stabilize a form of collinear spin order, before 
orbital order transitions simultaneously distort the cubic 
structure and cant spins at lower temperatures. C0V2O4 
is a unique case. Though a collinear ferrimagnetic tran¬ 
sition near Tni = 150 K^SHSSl } las i on g been observed 
and there have been recent reports suggesting a weaker 
spin canting transition near T^ 2 = 90K.221, the material is 
reported to retain cubic global symmetry to lowest mea¬ 
sured temperature} 27 } This last observation is especially 
interesting, as it contradicts all known models for this 
class of materials. We discuss it in more detail below. 

Enhanced spin-lattice coupling and associated phe¬ 
nomena resulting from orbital order have been reported 
on several occasions, and theory has predicted a number 
of different ordering patterns 3 ® 3 - 7 '. Given the difficulty 


of accessing the orbital order parameter experimentally, 
conclusions have been argued indirectly, and debate has 
centered on the relative importance of exchange, spin- 
orbit coupling and various crystal field terms in the mag¬ 
netic Hamiltonian. A secondary line of debate has fo¬ 
cussed on whether localized orbital models are appropri¬ 
ate at all, or whether one should approach the problem 
from the itinerant limit for some material^®!!. Indeed, 
the spinel vanadates were initially studied in the con¬ 
text of localized-itinerant crossover physics, and are pre¬ 
dicted to be metallic for V-V separations (Rw) l ess than 
the critical distance of R c = 2 . 97 AP Though known 
materials are all semiconductors, those with Rw ~ R c 
(e.g. MgV 2 0 4 , ZnV 2 04 and C0V2O4) have enhanced 
electron-lattice coupling and transport properties best 
described by the hopping of polaron} 33 } The activation 
energy of polarons is seen to decrease sharply (i.e. ma- 
terials become more metallic) with decreasing Rvv * 38 * 40 -* 
and increasing pressure 2 ^, but increase with disorder due 
to cation inversion 38 or mixed valence on the spinel B- 
sitp} C0V2O4, which has the shortest Rw, has been 
driven metallic over a limited temperature range with an 
applied pressure of 6 GPcP. 

Recent studies have expounded on the role played by 
the proximate itinerant state. One study showed that, 
although materials CdV204 and MnV 2 04 are well de¬ 
scribed by insulating models, the Mn compound exhibits 
anomalously high compressibility which has been a rgued 
to signify a breakdown of the Mott-Hubbard pictur^SEH. 
Materials with shorter Rw (A G {Fe, Zn, Mg, Co}) may 
lay in an intermediate phase, characterized by partial de- 
localization, mobility of large polarons and short-range 
magnetic ordei®. These arguments seem to be backed 
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by first principles calculations in C 0 V 2 O 4 , which sug¬ 
gest the coexistence of localized and itinerant electrons'® 
and electronic structure calculations in ZnV 2 04 , which 
explain the low-temperature transitions as a structural 
instability toward V-V dimerization, with no role for lo¬ 
cal orbital^H The role of orbital order in near-itinerant 
vanadates was further challenged by studies of the dop¬ 
ing series Mni_ x Co x V 204 , where x-ray and heat capacity 
measurements reveal the suppression of the low temper¬ 
ature tetragonal phase transition above critical doping 
x c ~ O.S® and neutron scattering reveals the decoupling 
of this transition from spin canting^- 

In light of the above discussion, the case of C 0 V 2 O 4 is 
particularly intriguing. This material is a direct analogue 
to MnV 2 04 , forming a near ideal spinel structure with 
V 3+ on the pyrochlore lattice and spin-only S=3/2 Co 2+ 
on the A-site. Its near critical V-V separation has made 
it the subject of several studies exploring t he role of the 
proximate itinerancy in 3d magnet systemd 27 -^MIlMll]. 
Though it is well-established that there exists a magnetic 
transition to a near colli near ferrimagnetic ordered state 
below T N1 ~ i 5 qk|M 129 | 38 | 431 , experimental situation 
at lower temperatures is quite unclear. While one study 
of C 0 V 2 O 4 powders reported a cusp in magnetization 
near T = 100 K and a sharp peak in heat capacity at T 
= 60 KpS, similar measurements of a single crystal sam¬ 
ple revealed only a magnetic anomaly at T = 75 K and 
no signature in heat capacity at any temperature below 
Tni ■ Neither study sa w a d eviation from cubic sym¬ 
metry down to T = 10 K^ 7 2Sl, which itself is surpri sing, 
as both orbital ordeiCSBEE and dimerization modeli®^® 
predict a symmetry lowering structural transition. Very 
recently, a neutron and strain study of crystals and pow¬ 
ders enriched with excess cobalt reported a spin canting 
transition Tn 2 = 90K, which they associated with the 
onset of an orbital glassiness, and a separate structural 
transition of unknown origin near T = 40 K 29 . 

The large variation in the literature is likely a reflec¬ 
tion of the combined weakness of observed effects and the 
enhanced sensitivity of this material to disorder. The 
differing heat capacity signatures in powders and crys¬ 
tals has been ascribed to the presence of random fields 
in cr ystalfP, as seen in FeCr 2 S 4 E3 , or could simply re¬ 
flect the difficulty in detecting first order transitions in 
heat capacity using the time-relaxation technique, as sug¬ 
gested elsewhere by uP®. Weak inversion, and in par¬ 
ticular cobalt on the spinel 13-site , is k nown to affect 
both the tem peratur e of transitions^ 9 ^ and the degree 
of localizatioipS122M2l i n this system. In fact, early magne¬ 
tization measurements noted a feature near 60 K, which 
seemed to disappear upon annealing and was attributed 
to the presence of C 02 VO 4 2 ®. Clearly, further measure¬ 
ments are needed before reaching final conclusions. 

In this article, we shed light on the outstanding ques¬ 
tions surrounding C 0 V 2 O 4 using a combination of neu¬ 
tron powder diffraction and inelastic scattering. Our 
diffraction results confirm the near cubic spinel structure 
with minimum cation inversion, and canted ferrimagnetic 


order below Tn = 156K. Significantly though, our anal¬ 
ysis further reveals the presence of a weak, first-order 
structural phase transition at T* = 90 K, the same tem¬ 
perature at which a spin canting transition is reported 
in Ref. 1291 This transition is evident in the tempera¬ 
ture evolution of both the lattice parameter, seen on long 
length-scales, and the oxygen atom positions, which has 
a short-range character. Our inelastic data show the suc¬ 
cess of local spin models in describing the magnon spec¬ 
trum, and suggest the presence of a small (1.25 meV) 
spin gap. Together, these observations strongly suggest 
that this material is most appropriately described by lo¬ 
cal orbital models, and indicate an orbital ordering tran¬ 
sition at T*. Proximate itinerancy is manifest primar¬ 
ily through a strong reduction of vanadium ordered mo¬ 
ments, canting angle and the magnitude of the structural 
distortion. 


II. MATERIALS PREPARATION AND 
EXPERIMENTAL DETAILS 

Powder samples of C 0 V 2 O 4 were prepared via solid 
state reaction at the National High Magnetic Field Lab in 
Tallahassee, FL, and extensive characterization data can 
be found on related samples in Refs. [27] and |43] As dis¬ 
cussed above, these references report a near cubic struc¬ 
ture to lowest temperatures and the formation of a net 
ordered moment below 150K, reflecting the ferrimagnetic 
spin order known from previous neutron measurement^®. 
Neutron scattering data in the current report was ob¬ 
tained entirely using instruments at Oak Ridge National 
Laboratory (ORNL) in Oak Ridge, TN. Powder diffrac¬ 
tion work was performed using the HB2A instrument at 
the High Flux Isotope Reactor (HFIRjfS®, followed by 
wide-angle time-of-flight (TOF) measurements using the 
Nanoscale-Ordered Materials Diffractometer (NOMAD) 
instrument at the Spallation Neutron Source (SNS) 51 . 
HB2A measurements were performed using ~3g of pow¬ 
der, placed in a vanadium can with helium exchange gas 
to ensure cooling. Patterns were measured independently 
with A = 1.538Aand A = 2.4lA neutrons from a germa¬ 
nium monochromator, and fit simultaneously. For NO¬ 
MAD measurements, about 100 mg of sample was mea¬ 
sured in 6 mm diameter vanadium cans for 20 minutes 
in an ILL Orange cryostat. NOMAD detectors were cal¬ 
ibrated using scattering from diamond powder, and the 
instrument parameter file for the Rietvelcl refinements 
was obtained from the measurements of the NIST stan¬ 
dard silicon powder. The structure factor S(Q) was ob¬ 
tained by normalizing diffraction data against a solid V 
rod, and subtracting the background using the IDL rou¬ 
tines developed for the NOMAD instrument. The pair- 
distribution function (PDF) was obtained by the Fourier 
transform of S(Q) with Q max = 31.5A -1 : 

G(r) = - [ Qma * Q(S(Q) - 1) sin (Qr)dQ (1) 

71 JQmin 
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Rietveld refinements were performed using Fullproi’ 5 -^, 
and PDF fits were performed using the PDFgui 
software®!. 

Inelastic neutron scattering (INS) measurements were 
performed using the hybrid spectrometer HYSPEC at 
the SNS. HYSPEC is a highly versatile direct geome¬ 
try spectrometer that combines TOF spectroscopy with 
the focusing Bragg optics 54 . The incident neutron beam 
is monochromated using a Fermi chopper with short, 
straight blades, and is then vertically focused by Bragg 
scattering onto the sample position by highly oriented 
pyrolytic graphite. HYSPEC employs 3He linear posi¬ 
tion sensitive tube detectors that are assembled into 20 
sets of 8 -packs that cover an angular range of 60° in the 
horizontal scattering plane and a vertical acceptance of 
15°. For the INS measurements, the C 0 V 2 O 4 powder 
was loaded in an aluminium can and placed in a closed 
cycle refrigerator capable of reaching a base temperature 
of 5 K. Data was collected with incident energies Ei = 
60 meV, 35 meV, 15 meV and 7.5 meV, and Fermi chop¬ 
per frequencies of 300 Hz, 180 Hz, 180 Hz and 420 Hz, 
respectively. 


III. NEUTRON SCATTERING RESULTS 
A. Powder Diffraction 

Neutron powder diffraction measurements were first 
performed using the HB2A instrument with two neutron 
wavelengths, and main results are shown in Figure |T| 
Figs, [l] (a.)-(c.) shows patterns from 1.54A measure¬ 
ments at T = 280 K, 115 K and 4 K, along with lines 
showing best-fit Rietveld refinements, performed using 
both wavelengths. Refinement of the pattern at T = 
280 K (Fig. [lja)) confirms the near ideal spinel struc¬ 
ture, with approximately 4% of the vanadium sites con¬ 
taining cobalt cations. Patterns at lower temperatures 
(Fig- 0b) -(c)) exhibit a visible increase in the intensity 
of several low-angle peaks associated with collinear ferri- 
magnetism, most prominently the cubic (111). Targeted 
measurement of this peak over a range of temperatures 
resulted in the data in Fig. [IJd), and simple power-law 
fits yield a transition temperature of Tn = 156.1 ± 0.5K 
and a power 2/3 = 0.68 ± 0.04. Consis tent with other 
neutron powder diffraction studieJ^®!, we saw no vis¬ 
ible intensity at the position of the cubic (200) Bragg 
peak associated with spin canting, yet refinements of the 
entire patterns to a model with Co 2+ spins along the 
c-axis and V 3+ spins allowed to cant along <110> from 
the antiparallel direction implied a small canting angle at 
base temperature. Specifically, at T = 4 K, best fits give 
ordered moments Mc 0 = 3.05(4)/r B , M v> || = 0.70(4)/r B , 
and My ± = 0.06(4)^ B , where M v j| and My,j_ denote or¬ 
dered vanadium moments antiparallel and perpendicular 
to Mco, respectively. The value for Mc 0 is within error 
equal to the expected full-moment value for S=3/2 Co 2+ . 
The fitted values for M V || and My ± within our model im- 
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FIG. 1: Powder diffraction patterns taken using the HB2A 
diffractometer with A = 1.54A at temperatures T = 280 K 
(a.), 115 K (b.) and 5 K(c.). Here, black dots are data and 
solid lines are the results of Rietveld refinements described in 
the main text. Arrows denote the position of the cubic (111) 
Bragg reflection, which serves as an order parameter for the 
collinear ferrimagnetic state. A power law fit to the temper¬ 
ature variation of this peak (d.) gives Tn = 156.1 ± OAK. 


ply a greatly reduced vanadium moment of My = 0.70/r B , 
with a canting of 5 ± 4°. Though the fitted canting angle 
is not significantly different from zero, this result is con¬ 
sistent with the reports of canting in Ref. [29] for Co-rich 
samples and Ref. 0S] for near itinerant members of the 
series Mn 2 _ x Co x V 204 . A complete listing of refinement 
parameters for each of the four measured temperatures 
is listed in Table [H 


B. TOF Diffraction and Pair-Distribution Function 
Analysis 

To search for local distortions of oxygen octahedra or 
other signatures of orbital glassiness, we supplemented 
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T (K) a (A) _ 

280 8.39976(8) 0.23965(7) 

115 8.39102(8) 0.23971(7) 2.43(4) 0.60(4) 0.04(4) 

85 8.39007(7) 0.23965 2.78(4) 0.67(4) 0.05(4) 

4 8.38924(7) 0.23964(7) 3.05(4) 0.70(4) 0.06(4) 


TABLE I: Parameters extracted from simultaneous Rietveld 
refinements of NPD patterns taken using A = 1.54A and 
2.4lA neutrons on the HB2A instrument at HFIR. Details 
about the fit model are provided in the main text. 


the above data with diffraction measurements using the 
NOMAD instrument, which has a significantly higher 
count rate and detector coverage more suitable for a PDF 
analysis. Data were analyzed using both Rietveld refine¬ 
ment in reciprocal space and with real-space PDF anal¬ 
ysis, described above, and the main results are shown 
in Figs. [2] and [3j respectively. Rietveld refinements used 
the same model and gave results largely consistent with 
conclusions from the HB2A analysis. Specifically, it was 
confirmed that the structure is very close to an ideal cu¬ 
bic spinel at all temperatures, and that spins order into 
the canted ferrimagnetic state described above. A rep¬ 
resentative pattern and associated fit line are shown in 
Fig.ga .). Additionally, the sizeable increase in signal-to- 
noise on the NOMAD instrument brings to light several 
interesting new features. Despite the success of the cubic 
spinel model across the full temperature range, Fig.pjb.) 
shows a clear, discontinuous jump of Aa = 0.0014W in 
the refined lattice parameter at a critical temperature T* 
= 90 K. The magnitude of this jump is within the error 
bars of previous x-ray measurements^, and this obser¬ 
vation is consistent with previous literature. Moreover, 
though the fractional position parameter of the oxygen in 
the Fd3m space group (Fig. |2](c.)) shows a smooth and 
continuous variation with temperature, the atomic dis¬ 
placement parameter (ADP), Bj so , which reflects the root 
mean square displacement of oxygen atoms from equilib¬ 
rium, shows a jump similar to lattice parameter at the 
exact same T*. Together, these observations strongly 
suggest the existence of an unresolved, first-order struc¬ 
tural phase transition at the limits of instrument resolu¬ 
tion. 

The existence of a structural transition at T* = 90 K 
is particularly important in light of previous polarized 
neutron scattering work on single crystals, which has 
established that a spin-canting transition exists at the 
exact same temperature 29 . Our own diffraction data is 
also not inconsistent with such a conclusion. In Table [TT| 
we show refined magnetic moments at selected tempera¬ 
tures assuming the canted state described above' 5 -^. One 
can see a smooth development of ferrimagnetism below 
Tjm, with moment sizes (Fig. [2](e.)) equal within error 
to HB2A fits and an implied canting angle for vanadium 
spins of 15 ± 5° at base temperature. This canting angle 
is statistically significant, though larger than the values 
in Table [T] by a factor of 3-4; this reflects the difficulty 
of reliably determining this quantity from powder refine- 


This canting angle is also a factor of 2 
larger than implied by single crystal measurements in 
Ref. [221 but notably equal to results from powder re¬ 
finements in the same study. In Fig. |2][f.), we show the 
temperature variation of My^, which is proportional to 
the cubic (200) Bragg peak intensity, used in Ref [29] to 
detect the spin canting transition in single crystals of 
this material. Though the magnitude of My ± (and thus 
canting angle) seems to grow with decreasing tempera¬ 
ture, statistical error bars in our data set are too large 
to discern a clear kink at T = T*. We do note, however, 
that the total increase of M?, over the measured tem- 
perature range is the same as seen in Ref. [291 and if one 
assumes a transition, fits to a power-law temperature de¬ 
pendence yield T* = 92K ± 5K. Thus, we conclude that 
our data is consistent with the previously reported spin 
canting transition at T*^. In light of this, the coincident 
structural transition seen in the current data set strongly 
suggests that this temperature can be associated with the 
ordering of V 3+ orbitals. 

The same data set was subject to a PDF analysis, us¬ 
ing the transform given by Eq. [TJ and assuming a cubic 
spinel structure. The main results, shown in Fig. [3j tell 
a similar story to the one above. The cubic spinel struc¬ 
ture describes the data well at all temperatures (see, e.g., 
Fig. [3](a.)), and there is a distinct jump in the inferred 
lattice parameter and the atomic displacement parame¬ 
ter of oxygen atoms at T* = 90K. Fits were performed 
over a number of increasingly shortened ranges in real- 
space: 1.5A< R R,„ ; ,x, with Rmax varying from 30A to 
15A. Additionally, one fit was performed using a range 
10 A< R < 45A, to explicitly suppress the effect of local 
correlations. As shown in Fig. [3](b.), the magnitude of the 
jump in lattice parameter at T* was found to be largely 
independent of fit range, consistent with expectations for 
a true long-range effect. In contrast, the displacement 
of the oxygens decreases monotonically with increasing 
Rmax and is strongly suppressed for the fit where local 
correlations are neglected. This suggests that the atomic 
displacement parameter of oxygens (Fig. [3](c.)) has an 
overwhelmingly local character. As this parameter rep¬ 
resents the root mean square positions of oxygens form- 


T (K) 

a (A) x 0 2 - 

Mco (mb 

) M v,\\ (Mb 

) M v ,± (mb ) 

144 

8.3909(2) 0.23968(5) 

2.02(4) 

0.53(5) 

0 

120 

8.3902(2) 0.23968(5) 

2.50(4) 

0.62(4) 

0.06(6) 

100 

8.3894(2) 0.23966(5) 

2.74(3) 

0.68(4) 

0.08(6) 

80 

8.3907(2) 0.23966(5) 

2.92(3) 

0.73(4) 

0.13(6) 

60 

8.3903(2) 0.23963(5) 

3.01(3) 

0.73(3) 

0.17(6) 

40 

8.3905(2) 0.23962(5) 

3.07(3) 

0.73(3) 

0.20(5) 

20 

8.3905(2) 0.23962(5) 

3.09(3) 

0.73(3) 

0.21(5) 

6 

8.3905(2) 0.23961(5) 

3.09(3) 

0.73(3) 

0.19(6) 


TABLE II: Fit parameters at selected temperature, extracted 
from Rietveld refinements of NPD patterns taken using the 
NOMAD instrument at the SNS. Details about the fit model 
are provided in the main text, and further parameters for the 
entire temperature range can be found in the Supplemental 
Materials. 


x 0 2- Mco (mb) M v n (mb) My,± {^b) ments alone. 
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FIG. 2: Main results from a Rietveld refinement of diffraction 
data taken with the NOMAD instrument. TOF data were 
fit to the ideal cubic structure and the canted spin model 
described in the main text. Panel (a.) shows a typical pattern 
from data collected using a single NOMAD detector bank. 
Here, the solid line shows the best fit to the above model. 
The temperature dependence shows a jump in the fit lattice 
parameter (panel (b.)) and atomic displacement parameter 
for oxygens (panel (c.)) at the same temperature, T* = 90 
K. Neither the fractional position parameter for the oxygens 
(panel (d.)), nor the ordered moment sizes (panel (e.)) show 
a signature at T*. Panel (f.) is the square of the transverse 
V 3+ ordered moment, equivalent to the cubic (200) Bragg 
peak, commonly thought to indicate spin canting. Solid line 
in this panel is a fit to a power-law temperature dependence, 
and implies T* = 92K ± 5K. 


ing the octahedra surrounding V 3+ octahedra, this ob¬ 
servation might suggest the onset of orbital glassiness, as 
previously conjectured^. As with the Rietveld analysis, 
we were able to fit the PDF data assuming a number of 
different structural distortions at low temperature, but 
none showed a significant improvement when compared 
to the cubic spinel structure. The weakness of the current 
effect does not allow us to comment further. 


C. Inelastic Neutron Scattering 


Further support for an orbital glass picture comes from 
inelastic neutron measurements on the HYSPEC spec¬ 
trometer, summarized in Figs. [4] and [5] Figs. [4j)a.) and 
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FIG. 3: Results of a PDF analysis of the NOMAD data, as¬ 
suming a cubic spinel structure and neglecting scattering from 
spin order. Panel (a.) shows a typical fit at base tempera¬ 
ture, demonstrating the ability of the model to describe the 
PDF patterns. Panels (b.) plot the refined lattice parame¬ 
ter as a function of temperature, as inferred from fits of real 
space data over several distinct data ranges. From top to 
bottom, data was fit using the range lA < R < R ma x, for 
R max = 15A, 20A, and 35A, whereas the bottom panel used 
10A < R < 45A. Panels (c.) show equivalent plots of the 
oxygen ADP. Both sets of panels reveal a discontinuity at 
T* = 90K, but exhibit a different dependence on fit range, as 
discussed in the main text. 


(b.) show the low-energy magnon spectra at T = 6 K, 
taken with Ej = 15 meV and 7.5 meV neutrons, respec¬ 
tively. Immediately apparent is the existence of a highly- 
dispersive, bright excitation mode emerging out of the 
position Q = 1.3A -1 , which is associated with the cu¬ 
bic (111) Bragg peak. Closer inspection reveals weaker 
excitations emerging from Q = l.sA^ 1 and Q = 2.lA _1 , 
the location of the (200) and (220) peaks, respectively. 
Figs. |4|c.) and (d.) confirm that these excitations give 
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way to a correlated paramagnetic phase at temperatures 
T > Tn- The flat band of scattering in Fig. [d|b.) is 
seen only with E; = 7.5 meV, but not with E; = 15 meV 
(Fig.Qa.)), 35 meV, or 60 meV (Fig.[5ja.)). For this rea¬ 
son, we conclude it is likely spurious, though its origin is 
unknown. 

Plots of scattering intensity versus energy transfer re¬ 
veal a small peak in scattering at A ~ 1.25meV above 
both magnetic Bragg peaks, which is suggestive of a spin 
gap. Though we are unable to confirm the complete ab¬ 
sence of scattering at lower energies due to experimental 
constraints, the magnitude of this gap is nearly a factor 
of 2 smaller than reported for Mn 2 _ x Co x V 204 for Co 2+ 
concentrations up to x = 0.8®. This is consistent with 
expectations for an orbital order state with strong renor¬ 
malization due to proximate itinerancy. 

In Fig. [5] we show a comparison of inelastic data taken 
with Ei = 60meV, and simulations based on a semiclas- 
sical spinwave Hamiltonian: 


n = J2 J bi s i • S J + Z D ' (Si ■ (2) 

<i,j> i 


along with the parameters in Table III Here,the sums 
are over all spins, including both cation sites, and in¬ 
teractions are truncated beyond nearest neighbors for 
the Co 2+ sites (nearest neighbor Co-V interactions) and 
next-nearest neighbors for the V 3+ sites (nearest neigh¬ 
bor V-V and nearest neighbor Co-V interactions only). 
We have used Jab to denote the exchange interaction be¬ 
tween A-site cobalt and B-site vanadium spins, while Jbb 
(J BB ) describes interactions between nearest-neighbor 
vanadium spins with the same (different) c-axis posi¬ 
tional coordinates. Single-ion anisotropy was set to zero 
for cobalt (Da ), consistent with the near-cubic structural 
symmetry, and fixed along local < 111 > directions for the 
vanadium sites (D B ). A similar analysis was successfully 
used by us to describe spin excitations in FeV 20 j 2 Uand 
in Ref. [151 to describe the Mni_ x Co x V 2 C >4 family of com¬ 
pounds. 

Notably, the material parameters listed in Table |III| 
and used to create Fig. §M were not extracted from fits 
to the inelastic data (Fig. |5](a.)), but rather were extrap¬ 
olated from parameters reported for single-crystalline 
MnQ 4 Coo.eV 2 O, and adjusted to account for Neel 
temperature, canting angle and spin gap. It should be 
pointed then out that the Heisenberg exchange parame¬ 
ters in Mno. 4 Co 0 . 6 V 20 4 are hugely anisotropic, as inter¬ 
actions between V-V pairs in that material are strongly 
modified by the existence of orbital order and the as¬ 
sociated structural distortion. Specifically, orbital order 
results in J BB and J BB that differ significantly not only 
in magnitude, but also sign. Though the lattice structure 
in C 0 V 2 O 4 is more isotropic, we found we were unable to 
satisfactorily describe the inelastic data in Fig. [5] without 
retaining the anisotropy in the Heisenberg exchange pa¬ 
rameters, and specifically the difference in sign between 
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FIG. 4: Powder inelastic scattering data, taken using the 
HYSPEC instrument at the SNS. Panels (a.) and (b.) show 
data taken at base temperature, T= 6 K, with Ei = 15meV 
and Ei = 7.5meV, respectively. The most obvious feature is 
the highly dispersive excitation, nearly touching the elastic 
line at the position of the (111) magnetic Bragg peak, and 
weaker bands are seen about (200) and (220). Vertical dotted 
lines denote the locations of the three aforementioned Bragg 
locations. Equivalent measurements taken at T = 200 K (pan¬ 
els (c.) and (d.)) show only correlated paramagnetic behav¬ 
ior, confirming that the above dispersive bands are magnetic 
in origin. Constant-Q cuts at the positions of the (111) and 
(200) (panels (e.) and (f.)) reveal a small bump, consistent 
with a spin gap of A = 1.25 meV in this material. 


J BB and J BB . The significance of this observation is dis¬ 
cussed in more detail below. 

Even without fitting, this basic spinwave model is able 
to account for several specific features of the C 0 V 2 O 4 
data, including the minima at 1.3 A -1 and 1.5 A” 1 , 
the ridge of intense scattering between 14 meV and 20 
meV, and the second ridge of scattering near 37 meV. Ar¬ 
guably, the largest failure of the model is an overestima¬ 
tion of scattering at higher Q, including the prediction of 
an intense mode near Q = 3-OA^ 1 . This observation indi¬ 
cates a deviation in C 0 V 2 O 4 from the assumed spherical 
approximation for the local magnetic form factor, and in 
fact may indicate that the atomic orbitals are distributed 
over larger distances in real space, as one might expect 
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FIG. 5: Plots of measured (left) and simulated (right) neutron 
inelastic scattering data with Ei = 60 meV, plotted on a log¬ 
arithmic scale. The simulated scattering data is the result of 
a spin-wave analysis, using parameters described in the main 
text and using spherical approximations for the local form 
factors of Co 2+ and V 3+ moments. Effects of instrument res¬ 
olution have not been included. 


for a material approaching an itinerant crossover. 


IV. DISCUSSION AND CONCLUSIONS 

Collectively, the neutron scattering results presented 
here draw a clear and consistent picture for the physics 
underlying C 0 V 2 O 4 . The most striking new result is the 
first order structural phase transition at T* = 90 K, which 
by itself addresses the largest outstanding question in the 
study of this material- namely the unexpected preserva¬ 
tion of the high symmetry cubic phase to lowest temper¬ 
atures. Here, we see that the predicted symmetry low¬ 
ering transition is not absent, but has merely been sup¬ 
pressed in magnitude to ^ ~ 10 -4 , presumably another 
effect of the near-itinerancy. The subtlety of the transi¬ 
tion and our inability to resolve details of the low sym¬ 
metry structure does not allow us to strictly determine 
whether it is driven by orbital ordei®^ or an instability 
towards homopolar bond formatior l — 1— 1 If one takes into 
account, however, previous reports of spin canting at the 

Ei Jab Jbb Jbb Da Db 
60 3 8 -4.59 0 -0.39 

TABLE III: Parameters used to produce simulated neutron 
scattering spectrum in Fig. [5])b) . All values are given in units 
of meV. 


same temperature 2 ® and the small spin gap seen in pow¬ 
der inelastic data at base temperature, the experimen¬ 
tal situation surrounding C0V2O4 is strongly reminis¬ 
cent of other ferrimagnetic spinels more unambiguously 
described by the Mott insulating picture 4 ®^!. Thus, it 
seems most appropriate to describe the 90 K transition 
in C0V2O4 as an ordering transition for localized V 3+ 
orbitals. In further support of this idea, we note that 
our successful description of inelastic scattering invoked 
a localized spin Hamiltonian, and used parameters which 
imply a strong anisotropy in Heseinberg exchange param¬ 
eters. This magnetic anisotropy echoes previous reports 
of an anisotropic compression in single crystali^. With¬ 
out a much larger tetragonal distortion than what we are 
reporting, the most natural way of explaining these ex¬ 
perimental signatures is though an emergent electronic 
anisotropy, such as one expects in an orbitally ordered 
state. 

On a final note, we point out that our data allows us 
to rule out any additional structural or spin transitions 
below T*, within the sensitivity of our probe. This is 
in direct contrast to reports of transitions at 40 I'PH, 60 
KpSl or 75 KP, seen in both powders and single crystals. 
We believe the variation on this point in the literature is 
largely a result of the difficulty in preparing pure samples 
of this material and its peculiar sensitivity to disorder. 
It is known that the most frequent impurity phase in the 
preparation of C0V2O4 is V2O3, which is strongly cor¬ 
related and has transition temperatures which are sensi¬ 
tive to level of defectsP. On the other hand, it is also 
known that the presence of excess cobalt on the vana¬ 
dium sit e stro ngly affects the degree of localization in this 
material 38 ] 39 ] 49 ^, drives changes in Neel temperature 3 ^, 
and has led to at least one report of a transition be¬ 
tween 45 K and 70 K which disappears upon annealing®. 
Moreover, it has been suggested that enhanced sensitivity 
to disorder as materials approach the itinerant crossover 
affects both the level of localization®! and might fur¬ 
ther explai n the m any reports of glassy correlations and 
dynamic' 2 ' —I— 1 — 1 . X-ray and neutron diffraction mea¬ 
surements on our sampies reveal no detectable impurity 
phases and a cation inversion (x ~ 0 . 04 ) that is among 
the lowest in the literature, which increases confidence in 
our results. A systematic multi-sample study will proba¬ 
bly be required, however, to disentangle complicated dis¬ 
order effects in this material and discern the true number 
of transitions. 

In summary, we have presented strong evidence for a 
structural phase transition in C0V2O4 at a temperature 
of T* = 90 K. This transition, coupled with previous re¬ 
ports of spin canting and the presence of a spin gap, lead 
us to believe that this material is best described by a 
local orbital ordering picture, though it is certainly pos¬ 
sible that ferrimagnetic bond formation is playing a role. 
We suggest that the main effects of the nearby itiner¬ 
ant state are the strong reduction of the magnitude of 
the structural distortion, the reduction of the vanadium 
ordered moment and canting angle associated with the 
















orbital order, and the modification of the magnetic form 
factor, indicative of a near delocalized spin distribution. 
Confirmation of these interpretations could come from 
future pressure experiments, where the known insulator- 
metal transition should completely suppress the reported 
structural transition and spin gap. 
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